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ABSTRACT 

Radial-velocity observations of Kepler candidates obtained with the SOPHIE and HARPS-N spectrographs have permitted unveiling 
the nature of the five giant planets Kepler-41b, Kepler-43b, Kepler-44b, Kepler-74b, and Kepler-75b, the massive companion Kepler- 
39b, and the brown dwarf KOI-205b. These companions were previously characterized with long-cadence (LC) Kepler data. Here we 
aim at refining the parameters of these transiting systems by i) modelling the published radial velocities (RV) and Kepler short-cadence 
(SC) data that provide a much better sampling of the transits, ii) performing new spectral analyses of the SOPHIE and ESPaDOnS 
spectra, after improving our procedure for selecting and co-adding the SOPHIE spectra of faint stars (^p > 14), and iii) improving 
stellar rotation periods hence stellar age estimates through gyrochronology, when possible, by using all the available LC data up to 
quarter Q17. Posterior distributions of the system parameters were derived with a differential evolution Markov chain Monte Carlo 
approach. Our main results are as follows: a) Kepler-4 lb is significantly larger and less dense than previously found because a lower 
orbital inclination is favoured by SC data. This also affects the determination of the geometric albedo that is lower than previously 
derived: Ag < 0.135; Z?) Kepler-44b is moderately smaller and denser than reported in the discovery paper, as a consequence of 
the slightly shorter transit duration found with SC data; c) good agreement was achieved with published Kepler-43, Kepler-75, and 
KOI-205 system parameters, although the host stars Kepler-75 and KOI-205 were found to be slightly richer in metals and hotter, 
respectively; d) the previously reported non-zero eccentricities of Kepler-39b and Kepler-74b might be spurious. If their orbits were 
circular, the two companions would be smaller and denser than in the eccentric case. The radius of Kepler-39b is still larger than 
predicted by theoretical isochrones. Its parent star is hotter and richer in metals than previously determined. 

Key words, planetary systems: individual (Kepler-39, Kepler-41, Kepler-43, Kepler-44, Kepler-74, Kepler-75, KOI-205) - stars: 
fundamental parameters - techniques: photometric - techniques: spectroscopic - techniques: radial velocities. 


1. Introduction 

Thanks to unprecedented photometric precision and temporal 
coverage, the Kepler space telescope has discovered over two 
thousand small-size d planetary candidates with radii 7?p < 4 /?0 
( [Burke et al.||2014[ ). At the same time, it has provided the ex¬ 
oplanet community with more than two hundred Jupiter-sized 
candidates, thus triggering further studies on the structure, for¬ 
mation, and evolution of giant companions as well as on their 
atmosphere, if the optical and/or nIR occultations are observed. 

Since 2010 we have been following up several Kepler giant 
candidates orbiting faint stars with Kepler magnitudes > 14 
using the SOPHIE spectrograph at the Observatoire de Haute 
Provence (Prance). In addition to determining the fraction of 
fal se pos itives among the Kepler giant candidates ( Santerne~et| 
[aL [ 2012 ] ), this intensive follow-up allo wed us to characteriz e 
the giant planets Kepler-4Ib/KOI-196b ( jSanterne et al.| 20 lT^, 


Kepler-43b/KOI-135b and Kepler-44b/KOI-204b (Bonomo et 
jaL] |20T^, Kepler-7 4b/KOI-200b and Kepler-75b/KOI-889b 
( jHebrard et al.|201^ ; the massive companion Kepler-39b/KOI- 
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423b in the brown-dwarf desert ( jBouchy et al. 201 1| ), which 
could be either an extremely massive plan et or a low-mass b rown 
dwarf; and the brown dwarf KOI-205b ( jPiaz et ar]|2013| ). Por 
two planets, that is for Kepler-74b and Kepler-75b, additional 
radial-velocity measurements were taken with the HARPS-N 
spectrograph ( jCosentino et al. 2012]), which has been installed 


at the Tel escopio Nazionale Ga 
2012 (see jHebrard et al.|2013| ). 


ileo at La Palma island in Spring 


All these giant companions were characterized using Kepler 
data with long-cadence (LC) temporal sampling of 29.42 min, 
usually because short-cadence (SC) photometric measurements, 
that is one point every 58 s, were not available at the moment of 
publication. However, the long-cadence sampling presents the 
strong inconvenience of distorting the transit shape. This effect 
leads to longer transit durations, more V-shaped transits, hence 
lower ratios between the semi-major axis and the stellar radius 
than the true ones. This yields lower stellar densities from the 
third Kepler law and thus makes both stellar and plan etary radii 
appear larger than t hey actually are ( [Kipping 2010[ ). To over¬ 
come this problem, [Kipping[ ( [2010[ ) suggested to perform the 
transit fitting by oversampling the transit model and then bin¬ 
ning the model samples to those of the LC before computing 
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the chi-square or the likelihood function. His Eq. (40) suggests a 
simple way to choose the resampling resolution, given the pho¬ 
tometric precision of the light curve. Following this prescrip¬ 
tion, [Kippmg^B^^ ( |2011| ) analysed the LC data of Kepler-4b 
through Kepler-8b using a bin number of 4. 

When modelling the transits of the Kepler planets observed 
with SOPHIE and HARPS-N, we therefore followed this sug¬ 
gestion by |Kippin^ ( |2010| ) and oversampled the transit model by 
a factor five, which is higher than recommended by his Eq. (40). 
However, in some cases, the analysis of short-cadence data is 
mandatory especially when the orbital period is close to an inte¬ 
ger multiple of the LC sampling because this prevents the 
transit from being well sampled in orbital phase. The most evi¬ 
dent case is K epler-43b, whose orbital period is P = 147.99 -^Pic 
(see Fig. 5 in |Bonomo et al. [2012a| ). In addition, the massive 
companion Kepler-39b did not present an optimal coverage of 
the transit i ngress and egress in th e quarters Q1 and Q2 that were 
analyse dby|Bouchy et al.|(|2011| ) because P = 1032.14-^ric (see 
Fig. 9 in |Bouchy et al.|201^ 

Moreover, by performing an homogeneous analysis of tran¬ 
sit photometry from space and oversampling the Kepler LC data 
by a factor of ten instead of five, |Southworth| ( |2012| ) derived 
smaller stellar and planetary radii for some of the aforemen¬ 
tioned giant companions, although his results agree with ours 
within 2 cr. This would indicate that an oversampling of the tran¬ 
sit model by a factor of five might not be ideal in all cases, hence 
an independent analysis of SC data is certainly recommended. 
Nevertheless, some of the slightly different results obtained by 
|Southworfii| ( |2012| ) are also due to a better ephemeris and tran¬ 
sit signal-to-noise ratio (S/N) because, in most cases, he used 
longer temporal series than we did, up to quarter Q6. 

In this paper, we report the results of our analysis of the 
Kepler SC data of Kepler-39, Kepler-41, Kepler-43, Kepler-44, 
Kepler-74, Kepler-75, and KOI-205, along with the previously 
published radial velocities (RV). Kepler LC data up to quar¬ 
ter Q17 were used to refine stellar rotation perio ds by means 
of both generalised Lomb-Scargle periodograms ( |Zechmeister| 
|& Kurster|2009| ) and autocorrelation functions, when an unam¬ 
biguous peak with FAP < 0.01% could be identified. This allows 
us to estimate system ages through gyrochronology ( [Mamajek] 
|& Hillenbrand||20()^ , after deri ving the B-V index colour and 
its uncertainty from Eq. (3) in [Sekiguchi & Fukugita| ( |2000| ). 
The IDs and Kepler magnitudes of the parent stars are listed in 
Table This work aims at refining the characterization of these 
systems and possibly clarifying the apparently unusual proper¬ 
ties of Kepler-39b and Kepler-4 lb. Indeed, the former was found 
to have a larger r adius than predicted by theoretic al isochrones 
of |Baraffe et al.| ( |2003) ), and |Bouchy et al.| ( [2011| ) were unable 
to find any reasonable explanation for this behaviour. The latter 
seemed to be a non-infiated planet despite its proximity to the 
host star, and to occupy an atypical position in the radius-mass 
and radius-Teg diagra ms of giant planets (see Figs. 9 and 10 in 
[Santerne et al.|20lTa| ). 

Moreover, we performed new spectral analyses of the planet¬ 
hosting stars after improving our procedure of selecting, treat¬ 
ing, and co-adding the SOPHIE spectra. A revision of the at¬ 
mospheric parameters may have a significant impact on stellar, 
hence planetary, parameters. 

We recognize the merit of an approach to revisiting stellar 
or planetary parameters o f transiting systems that encompasses 
much larger samples (e.g., Torres et al.|2012t['Southworth|2012| ) 
than the one presented here. Our work differs in that it performs a 
self-consistent re-analysis taking into account both photometric 
and spectroscopic measurements and constraints within a coher¬ 


ent Bayesian framework to derive the posterior density distribu¬ 
tions of the full set of system parameters. 


2. Data 

2.1. Kepler photometry 

Short-cadence measure ments obtained with t he simple-aperture- 
photometry pipelin^ ([Jen kins et al. 2010 1 were downloaded 
from the MAST archiv^ Eleven quarters of SC data (Q3- 
Q7 and Q10-Q15) are available for Kepler-43; six quarters 
(Q10-Q15) for Kepler-75; four quarters (Q4-Q7) for Kepler-41, 
Kepler-44, Kepler-74, and KOI-205; and three quarters (Q12- 
Q14) for Kepler-39. 

The medians of the errors of SC measurements of Kepler- 
39, Kepler-41, Kepler-43, Kepler-44, Kepler-74, Kepler-75 and 
KOI-205 are 1.20 • 10"^ 1.27 • 10“^ 8.9 • lO'^, 2.02 • 10“^ 
1.24 • 10“^, 3.26 • 10“^, and 1.28 • 10“^ in units of relative flux, 
respectively. 

For all targets, the flux excess originating from background 
stars that are located within the Kepler photometric mask was 
subtracted separately for each quarter by using the estimates pro¬ 
vided by the Kepler teanj^ Indeed, this contamination of the tar¬ 
get flux dilutes the transits, making them appear shallower than 
they are, even though it usually does not exceed 5-7% of the total 
collected flux. 

All the SC data were used to model the transits of the seven 
giant companions. The signal-to-noise ratios of the phase-folded 
transits are of at least - 450 and, in some cases, higher than 
1 000. Thanks to these high S/N, we were able to derive stel¬ 
lar and planetary radii with uncertainties < 3% almost in all 
cases (see Tables [2][^. At this level of precision, errors on plan- 
etary radii are dominated by the uncertainties on stellar models 
( South worth|20 11 1 [20 1 2| ) and/or on the orbital eccentricity e and 
argument of periastron oj. Indeed, the uncertainties on e and oj 
from RV observations propagate into the transit parameter ajR^ 
and thus into the stellar density from the third Kepler law. Stellar 
density is then used as a proxy for luminosity to determine stellar 
hence planetary parameters (e.g., |Sozzetti et al.|2007| ) when no 
constraints from asteroseismology are available. This means that 
in our particular cases the additional use of LC data practically 
does not yield any significant improvement on planetary param¬ 
eters while introducing possible covariances between transit pa¬ 
rameters ( [Price & Rogers ||2014| ). For this reason, as mentioned 
before, we used LC data only to derive stellar rotation periods. 


2.2. Radial-velocity data 

The RV observations considered in this work are those listed 
in the announcement papers because no additional observa¬ 
tions with either SOPHIE or HARPS-N were carried out for 
these targets. The SOPHIE measurements were performed in 
high-efficiency mode with a resolution of - 40 000, and expo¬ 
sure times not exceeding 1 hr. The observations of KOI-200 
and KOI-889 carried out with HARPS-N were taken in high- 
resolution mode (the only available one) with a resolving power 
of - 110000 and exposures of 45 min and shorter than 25 min, 
respectively. Both SOPHIE and HARPS-N measurements were 


^ http://keplergo. arc.nasa.gov/PyKEprimerLCs. shtmlp 
^ http://archive.stsci.edu/kepler/data_search/search.php 
^ http://archive.stsci.edu/kepler/kepler_fov/search.php; fo r Season 2 
data of KOI-205, the value of crowdin g factor as derived by |Dfaz et al.| 


( 2013| ) was used (see Diaz et al. 2013| l. 
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Table 1. IDs, coordinates, and magnitudes of the planet-hosting stars Kepler-39, Kepler-41, Kepler-43, Kepler-44, Kepler-74, 
Kepler-75, and KOI-205 


Kepler name 

Object 

Kepler ID 

2MASS ID 

Kepler magnitude Kp 

Kepler-3 9 
KOI-423 
9478990 

19475046+4602034 

14.33 

Kepler-41 

KOI-196 

9410930 

19380317+4558539 

14.46 

Kepler-43 

KOI-135 

9818381 

19005780+4640057 

13.96 

Kepler-44 

KOI-204 

9305831 

20002456+4545437 

14.68 


Kepler name 

Object 

Kepler ID 

2MASS ID 

Kepler magnitude Kp 

Kepler-74 

KOI-200 

6046540 

19322220+4121198 

14.41 

Kepler-75 

KOI-889 

757450 

19243302+3634385 

15.26 

KOI-205 

7046804 

19415919+4232163 

14.52 



performed in obj_AB observing mode with fibre A centred on 
the target and fibre B on the sky. When needed, the obs ervations 
were corrected for moonlight pollution, as described inlBonomo 
|etaLl ( [2QT0| . 

2.3. Spectra 

The atmospheric parameters of the host stars, along with the stel¬ 
lar density derived from the transit fitting, are of fundamental 
importance for deter mining stellar, hence planetary, parameters 
( [Sozzetti et al.|2QQ7| ). 

While radial-velocity measurements can accommodate low 
S/N spectra, spectral analysis is more challenging. Indeed, some 
diffuse light in the SOPHIE spectrograph might affect the spec¬ 
tra at very low S/N acquired in high-efficiency mode. For this 
reason, we recently improved our procedure of treating and se¬ 
lecting the SOPHIE spectra to determine stellar atmospheric pa¬ 
rameters. In particular, spectra with an S/N lower than 14 were 
excluded from the co-addition. Those acquired in the presence 
of the Moon were corrected for the moonlight contamination by 
subtracting the background as estimated from fibre B. As is usu¬ 
ally done, the individual exposures were then set in the rest frame 
and co-added in a single master spectrum. 

In the case of Kepler-39, the co-added spectrum obtained this 
way shows deeper lines than simply co-adding all the SOPHIE 
spectra, as was previously done by |Bouchy et al.| ( [20111 ). This 
has a significant impact on the derivation of the atmospheric pa¬ 
rameters (see Sect. |4.1| ). 

The S/N of the SOPHIE master spectra at 600 nm and per 
element of resolution ranges between 100 and 170 for Kepler- 
39, Kepler-41, Kepler-43, Kepler-44, and Kepler-74. It is equal 
to 72 and 65 for KOI-205 and Kepler-75, respectively. 

Two host stars, namely KOI-205 and Kepler-39, were also 
observed with ESPaDOnS at the 3.6-m Canada-France-Hawaii 
Telescope in Mauna Kea as part of a programme dedicated to 
the characterization of Kepler planet-hosting star^ The objec¬ 
tive is indeed to carry out a better spectral analysis of the parent 
stars with a spectrograph that offers both a higher spectral res¬ 
olution (R ^ 65 000) and an extended spectral coverage (370 
- 1000 nm). These two targets were observed in ’object-fsky’ 
mode. The ESPaDOnS spectrum of KQI- 205 with a 5 /A ~ 90 
was previously used by |Diaz et al.| ( 2013 ) to determine the host 
star and brown dwarf parameters. Kepler-39 was observed with 
ESPaDOnS on September 28, 2012 and December 1, 2012, in 
a series of five exposures of ~ 40 min. The individual spec¬ 
tra as reduced by the CFHT Upena/Libre-Esprit pipeline were 
co-added after they were set in the rest frame and resulted in a 


master spectrum with S/N of 65 in the continuum at 600 nm per 
resolution element. This co-added spectrum is analysed here for 
the first time. 


3. Data analysis 

3.1. Spectral analysis 

To determine the effective temperature (Teff), surface gravity 
(logg), and iron abundance [Fe/H], the co-added SOPHIE spec¬ 
tra obtained with our improved selection and those acquired with 
ESPaDOnS were re analysed following the same procedures de¬ 
scribed in detail by |Sozzetti et al.| ( |2004[ |2006| ) and references 
therein. A set of ~ 60 relatively weak lines of Fe I and 10 of 
Fe II were selected, and EWs were measured using the TAME 
software ( [Kang & Lee|2012| ). Metal abundances were derived as¬ 
suming local thermodynamic equilibrium (LTE ), using the 20 10 
version of the spectral synthesis code MOOG ( |Sneden|1973| ), a 
grid of Kurucz ATLAS plane-parallel model stellar atmospheres 
( Kurucz||1993| ), and imposing excitation and ionisation equilib¬ 
rium. Uncertainties in the parameters were estimated following 
the prescriptions of jNeuforge & Magain| ( |1997) and [Gonzalez & 
Vanture[ ( ri998[ ) and rounded to 25 K in Tgff and 0.05 dex in log g. 

The derived atmospheric parameters are compared in 
Sect. [4.1 [ with the literature values, which were obtained with i) 
the iterative spectral synthesis package VWA ( [Bruntt et al.[2010| ) 
for Kepler-39, Kepler-41, Kepler-74, Kepler-75, and KOI-205, 
and u) the 2002 version of the MOOG code with the methodol¬ 
ogy described in [Bonomo et al.[ ( [2012a[ ) and Mortier et al. ( [2013 ) 
for Kepler-43 and Kepler-44. 


programme 12BF24, PI: M. Deleuil 


3.2. Combined analysis of Kepler and radial-velocity data 

To derive system parameters, a Bayesian analysis of Kepler SC 
photometry and radial-velocity measurements was performed, 
using a differential evolution Markov chain Monte Carlo (DE- 
MCMC) method ( [Ter Braak[2006[[Eastman et al.[2013[ ). For this 
purpose, the epochs of the SOPHIE and HARPS-N observations 
were converted from B JDutc into B JDtdb ( [Eastman et al.[20I(}| ), 
which is the time stamp of Kepler data. 

The transit fitting was performed using the model of 
Gimenez ( [2006[ [2009[ ). For this purpose, each transit was nor¬ 
malised by locally fitting a slope to the light-curve intervals of 
twice the transit duration before its ingress and after its egress. 
For Kepler-39, a linear function of time did not provide a satis¬ 
factory normalisation because of the short-term stellar variabil¬ 
ity (see Sect. [4.2[ ), hence a quadratic func tion of time was u sed. 
Correlated noise was estimated following Pont et al. ( 2006[ ) and 
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Bonomo et al. ( 2012b| ), and added in quadrature to the formal er¬ 


ror bars. However, it turned out to be very low, generally lower 
than one fifth of the formal photometric errors, as expected for 
high-precision space-b ased photometry (e.g., Aigrain et al.|2009 
[Bonomo et al.|201^ . 

Our global model has 12 free parameters when i) an ec¬ 
centric model was considered and ii) RV were taken with only 
one instrument (SOPHIE): the transit epoch To; the orbital pe¬ 
riod P\ the systemic radial velocity V/, the radi al-velocity semi¬ 
amplitude K; ^Je costu and sinoj (e.g., [Anderson et al. 


2011|); an additive RV jitter term Sj to account for possible jitter 
in the RV measurements regardless of its origin, such as instru¬ 
mental effects, stellar activity, additional companions, etc.; the 
transit duration from first to fourth contact Tu; the ratio of the 
planetary-to-stellar radii Rp//?=,; the inclination i between the or¬ 
bital plane and the plane of the sky; and the two limb-darkening 
coefficients (LD C) q\ = {Ua + ui)f and q 2 = 0.5ua/(Ua + ut) 
( |Kipping|2013| ), where u. and u\^ are the coefficients of the limb- 
darkening quadratic la^^jJ Two additional parameters, that is the 
HARPS-N systemic radial velocity and jitter term, were fitted 
when HARPS-N data were obtained as well (Kepler-74 and 
Kepler-75). Uniform priors were set on all parameters, in par¬ 
ticular with bounds of [0, 1] for q\ and q 2 ( Kipping [2013] ), lower 
limit of zero for K and s -^, and upper bound of 1 for e (the lower 
limit of 0 simply comes from the choice of fitting cos tu and 
^fe sintu). 

The posterior distributions of our free parameters were de¬ 
termined by means of our DE-MCMC code by maximising a 
Gaussian likelihood (see, e.g., Eq. 9 and 10 in Gregory ||2005| ). 
Eor each target, a number of chains equal to twice the number 
of free parameters were run simultaneously after being started at 
different positions in the parameter space but reasonably close to 
the system values known in the literature and/or obtained with an 
independent fit that was previously performed with AMOEBA 
( jNelder & Meadp965| ). The jumps for a current chain in the 
parameter space were determined f rom the other ch ains, ac¬ 
cording to the prescriptions given by [Ter Braak ( 2Q06| ), and the 
Metropolis-Hastings algorithm was used to accept or reject a 
proposed step for each chain. 

Eor the convergence of the chains, we required the Gelman- 
Rubin statistics, R, to be lower than 1.03 for all the parameters 
( Gelman et al. 12003] ) . Steps belonging to th e bum-in phase were 
identified following [Knutson et al. ( 2009[ ) and were excluded. 
The medians of the posterior distributions of the fitted and de¬ 
rived parameters and their 34.13% intervals are reported as the 
final values and their Icr error bars. When the distributions of 
the eccentricity and the RV jitter were found to peak at zero, we 
provided only the Icr upper limits estimated as the 68.27% con¬ 
fidence intervals starting from zero. Indeed, the medians of these 
distributions might yield misleading non-zero values. 

Einally, the Yonsei-Yale evolutionary tracks ( [Demarque et[ 
ar][2004[ ) for the effective temperature, metallicity, and density 


of the host stars were used to determine the stellar, hence com¬ 
panion, parameters ( [Sozzetti et al.|2QQ7} [Torres et al.|2Q12] ). 

Eor Kepler-41 and Kepler-74 only a circular model was 
adopted for the following reasons: the RV curve of Kepler-41 
is slightly asymmetric with respect to a sinusoid, very likely be¬ 
cause of residual effects from the correction of moonlight con¬ 
tamination and/or low S/N (from 13 to 20) spectra. These ar¬ 
tificial asymmetries in the RV curve tend to bias the solution 


^ = 1 - Wa(l - jd) - Wb(l - where 7(1) is the specific 

intensity at the centre of the disc and jd - cosy, y being the angle 
between the surface normal and the line of sight. 


of system parameters towards a low but non-zero eccentricity 
e ~ 0.14 even when including in the global fit the secondary 
eclipse that indicates that e cos oj is consis tent with zero (e.g., 
[Santerne et al.|201 la} [Quintana et al.|2013[ ). This would reduce 
the orbital configurations of a possible eccentric orbit to tu = 90 
or 270 deg. However, the expected circularization timescale is 
shorter than 100 Myr by assuming a modified tidal quality fac¬ 
tor of 2p = 10^ for the planet because of its short orbital period 
P = 1.85 days (hence small semi-major axis a = 0.031 au) and 
relatively low mass Mp ~ 0.6 Mjup for a Jupiter-sized planet. 
Indeed, there are no planets with mass comparable to Kepler-41 
and R < 3 d with a significant eccentricity. Eor these reasons, we 
adopted a circular model for Kepler-41. 

Eor Kepler-74, our DE-MCMC chains did not converge to¬ 
wards a unique solution when we varied the eccentricity, which 
resulted in very low acceptance rates. This was also noticed by 
[Hebrard et al.[ ( [201^ , who imposed a Gaussian prior on the or¬ 
bital eccentricity, solely based on a RV fit, in their combined 
analysis of Kepler and RV data. However, this prior inevitably 
affects the posterior distributions of orbit and transit parameters. 
Instead, we preferred to use a circular model given that current 
RV data evidently do not allow us to constrain the orbital eccen¬ 
tricity well. 

The Kepler-39 system parameters were obtained by using 
both eccentric and circular models because the 2(t significance 
of the eccentricity e = 0.112 + 0.057 cannot exclude that the 
orbit is perfectly circular ( [Lucy & Sweeney|1971[ ). 


4. Results 


4.1. Stellar atmospheric parameters 

The Teff and [Ee/H] of Kepler-41, Kepler-43, Kepler-44, and 
Kepler-74, which were determined with the procedure described 
in Sect. [3.1[ are consistent within 1 cr with the literature values. 
The surface gravities of Kepler-43 and Kepler-44 were found to 
be log g = 4.4 + 0.10 (Kepler-43) and log g = 4.1 +0.10 (Kepler- 
44), wh ich are lower than previously found by [Bonomo et al.| 
( [2012a[ ), that is logg = 4.64 ± 0.103 and 4.59 ± 0.14, respec¬ 
tively. These newly determined values are more consistent with 
the photometrically derived log g (see Tables]^ and [^. 

The st ar KQI-205 was f ound to be slightly hotter than re¬ 
ported by Diaz et al. ( [2013[ ), with Teff = 5400 ± 75 K, metal¬ 
licity [Ee/H]=0.18 ± 0.12, and logg = 4.7 ± 0.10, both with the 
ESPaDOnS and the co-added SOPHIE spectrum (see Table [^. 
The slightly hotter temperature has only a minor, almost negli¬ 
gible, influence on system parameters. 

The most striking differences with previously determined at¬ 
mospheric parameters were found for Kepler-39 and Kepler-75. 
The former is significantly richer in metals and slightly hotter 
than reported by [Bouchy et al.[p011[ ): Tgff = 6350 ± 100 K, 
[Ee/H] = 0.10+0.14, and logg = 4.4+0.15, to be compared with 
the previous estimates Teff = 6260+140, [Ee/H] = -0.29 + 0.10, 
and logg = 4.1 + 0.2. Very consistent values were derived from 
the analysis of the ESPaDOnS spectrum. This difference comes 
from the better selection and treatment of the SOPHIE spectra 
before co-adding the individual spectra (see Sect. [2.3[ ). It is not 
due to the different sp ectral analysis technique that was used by 
Bouchy et al.[ ( [2011] ), that is, the VWA package ( Bruntt et al. 
2010[ ). Indeed, when run on the new co-added SOPHIE spec¬ 
trum, VWA provided results that are almost identical to those ob¬ 
tained with MOOG: = 6360 + 100 K, [Ee/H] = 0.12 + 0.14, 

and log g = 4.4+ 0.3. These new values have an impact on stellar 
mass and age (see Table |^. 
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Kepler-75 is also found to be richer in metals than previously 
thought. Its atmospheric parameters determined with MOOG are 
Jeff = 5200+100 K, [Fe/H] = 0.30+0.12, and log g = 4.6+0.15. 
While the is consistent within 1 cr with the value derived by 
Hebrard et al.| ( 201^ , the difference in metallicity is about 2 a 


(see Table[7]). However, at these temperatures, the system param¬ 
eters of Kepler-75 do not significantly change as a consequence 
of the higher metallicity. 


4.2. System parameters 
4.2.1. Kepler-39 


Orbital and transit parameters obtained with our Bayesian DE- 
MCMC analysis in the eccentri c case agree well with those that 
were previously determined by |Bouchy et al. ( |2011| ). However, 
we found a lower significance of approximately 2 cr for the or- 
bital e ccentricity, that is ^ = 0.112 + 0.057, while [Bouchy et al.] 
( 2011| ) reported e = 0.122 + 0.023. Our larger uncertainty on 
the eccentricity indicates that it might be spurious, according to 
[Lucy & Sweeney| ( |l971 ). The Bayes factor of ~ 3 between the 


eccentric and the circular model, which was computed by using 
the truncated posterior mixture method ( |Tuomi & JonespoT^ , 
does not provide s trong enough evidence either for an eccentric 
orbit, according to |Kass & Raftery| ( |l995| ). 

As discussed in Sect. |4.1[ the stellar atmospheric parameters 
were refined thanks to both a better treatment of the SOPHIE 
spectra and a new ESPaDOnS spectrum. Specifically, a moder¬ 
ately hotter temperature Teff = 6350 + 100 K and a signific antly 
higher metallicity of 0.10 + 0.14 were found (see Sect. |4.1| ). 
For these new atmospheric parameters and the transit density, 
the Yonsei-Yale evolutionary tracks indicate a more massive and 
younger star: M* = 1.29 ^qq 7 M©, = 1.40 + 0.10 R©, and 

age of 2.1 ^q^ Gyr. The corresponding mass, radius, and density 
of Kepler-39b are Mb = 20.1+};^ Mj^p, = 1.24+^09 Rj^p, and 
Pb = 13.0^2 2 S cm“^. The se companion parame ters are consis¬ 
tent with those reported by |Bouchy et al.|p011| ), except for the 
stellar age, which is about half the value found by these authors 
(see Table 1^. Interestingly, this updated value of the age agrees 
well w ith the gyrochronology estimate ( [Mamajek & Hillenbrand| 


2008 

7 ^ 


that is t. 


gyr 


_ QJ+0.9 


Gyr, for the stellar rotation period 
= 4.50 + 0.07 d inferred from the Kepler LC light curve. 

Figure shows the position of Kepler-39b (empty red cir¬ 
cle) in the radius-mass diagram of transiting companions with 
masses between 10 and 100 Mjup for th e eccentric case. The 
dashed lines from top to bottom show the Baraffe et al. (|2003|) 
isochrones for 0.5, 1, 5, and 10 Gyr. As already noted by |Bouchy| 
|et al.| ( [MTT] ), the companion radius given by the eccentric solu¬ 
tion would be incompatible with theoretical isochrones with a 
probability > 95% (2 cr). However, the solution these authors 
proposed to explain the large radius of Kepler-39b, that is an in¬ 
creased opacity in the companion atmosphere, which previously 
seemed unlikely for the low stellar metallicity, now might apply 
for Kepler-39b. Indeed, our new estimate of [Fe/H] is signifi¬ 
cantly higher. 


As previously discussed, the eccentricity of Kepler-39b 
might be spurious according to the Lucy-Sweeney criterion. In 
the circular case, stellar and companion radii are slightly smaller 
than in the eccentric case, which implies a higher bulk density 
of 17.4^} 4 g cm“^ (see Table [^. The system age would be 
1.0^Q2 Gyr. The position of Kepler-39b in the circular case is 
shown in Fig. with a red filled circle. Figure displays the 



Fig. 1. Radius-mass diagram including transiting companions 
more massive than 10 Mjup. The dashed curves are the [Baraffe 


et al.| ( |2003 ) isochrones for 0.5, 1, 5, and 10 Gyr (from top to 


bottom). Both circular (filled red circle) and eccentric (empty 
red circle) solutions are shown for Kepler-39b. Blue labels indi¬ 
cate the brown dwarfs that were charac terized thanks to SOPHIE 
spectroscopic measurements (see aIso|Moutou et aI.|2013l|Diaz 
let al.|2014| ). 


phase-folded transit and RV curve and the best solutions ob¬ 
tained with both eccentric and circular models. 


4.2.2. Kepler-41 

The new system parameters deter mined with SC data si gnifi- 
cantly differ from those obtained by jSanterne et al. ( 201 la| ). The 
most important difference in the fitted parameters is found for the 
orbital inclination that, in turn, affects the determination of alRj, 


and the stellar density. Indeed, while jSanterne et al.| ( |2011a| ) 
found i = 88.3 + 0.7 deg, our new solution points to a consid¬ 
erably higher impact parameter with i = 82.51 + 0.09 deg (see 
Table[^. _ _ 

By analysing only LC data, |Southworth| ( |2012| ) found two 
solutions for the orbital inclination, one with i ~ 90 deg and the 
other one with i ~ 80-82 deg. He opted for the first because “the 
/ ~ 80 - 82 deg family occurs mainly for LD-fixed light-curve 
solutions, and results in weird physical properties”. However, by 
letting the LDC vary, our DE-MCMC run on SC data always 
converged toward the latter value of i. This occurred even when 
the chains were all started at values close to i ~ 90 deg. This 
ambivalence shows that, in some cases, fitting the LDC with LC 
data may lead to local minima that do not represent the true so¬ 
lution. 

Long-cad ence measurements were also used by |Quintana 
et al.| ( |2013| ) to validate this planet by analysing the phase 
curve. Curiously, these authors found an orbital inclination of 
85.4 ^q^ deg, which is in between the other two solutions. 
However, we point out that |Quintana et al.] ( 2013| ) were more in¬ 
terested in analysing the phase curve and did not explicitly men¬ 
tion which oversampling they adopted to model the LC transits. 

Our solution with low orbital inclination i = 82.5 deg is 
physically acceptable because the lowest possible inclination for 
the transit of Kepler-4 lb is ~ 78 deg. The corresponding stel¬ 
lar density ind icates a star that is lar ger and older than previ¬ 
ously found by jSanterne et al.| ( [2011a ) for the same atmospheric 
parameters. The derived logg = 4.278 + 0.005 is now more 
consistent with the spectroscopic value logg = 4.2 + 0.10 (see 
Table [^. A larger stellar radius implies a larger planetary radius 
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Fig. 2. Left panel. Top: phase-folded transit of Kepler-39b along with the transit models for the circular (red 
solid line) and eccentric (red dotted line) orbits. The two models are indistinguishable. Middle: residuals 
of the circular orbit. Bottom: residuals of the eccentric orbit. Right panel. Top: phase-folded radial-velocity 
curve of Kepler-39 and, superimposed, the Keplerian models for the circular (red solid line) and eccentric 
(red dotted line) orbits. Middle: 0-C of the circular orbit. Bottom: 0-C of the eccentric orbit. 


/?p = 1.29 + 0.02 R: 


jup, which now makes this planet more simi¬ 
lar to the other close-in hot Jupiters. Indeed, fro m the planetary 
parameters that were previously determined by [Santerne et al 


( 2011a| ) and Southworth| ( 2012| ), this object appeared quite rare, 
meaning: non-in flated despite its vicin ity to the parent star (see 
Figs. 9 and 10 in |Santerne et al.|2011a| ). 

Figure [^displays the phase-folded transit and radial-velocity 
measurements along with the best solution. 

Last but not least, the lower inclination also has an impact 
on ajR^ that becomes equal to 50.26 + 0.36 and, in turn, affects 
the determination of the geometric albedo (cf., e.g., Eq. (14) in 


Rowe et al. 2006 ). The latter was found to be considerably higher 


than the majority of ho t Jupiter s with measured optical occulta- 
tions by|Santerne et al.|(|2011a]), that is Ag = 0.30 + 0.07. On the 
contrary, our new solution with SC data indicates a significantly 
lower geometric albedo of Ag < 0.135 from the most recent 
value of the secondary-eclipse depth ( [Angerhausen et al.|2014| ), 
in agreement with theoretical expectation s for atmospheres o f 
hot Jupiters without scattering clouds (e.g.. Burrows et al. 2008| ). 
Figure]^ shows the albedo values and their Icr uncertainties as a 
function of the planet day-side equilibrium temperature Teq. The 
vertical dashed lines indicate the values of Tgq assuming perfect 
heat redistribution (left) or no redistribution in the atmosphere 
(right). 



1600 1800 2000 2200 2400 2600 

Day-side equilibrium temperature [K] 


Fig. 4. Geometric albedo of Kepler-4 lb as a function of its day- 
side equilibrium temperature. Tgq. The vertical dashed lines indi¬ 
cate the values of Teq assuming perfect heat redistribution (left) 
and no redistribution in the atmosphere (right). The grey band 
shows the albedo values allowe d by the Icr uncertaint y on the 
occultation depth determined by [Santerne et al. ( |2011a| ). 


4.2.3. Kepler-43 

System parameters derived from our DE-MCMC analysis are 
listed in Table and Fig. shows the phase-folded SC transit 
and radial-velocity curve along with the best-fit model. System 
parameters generally agree within 1 a with those determined by 
[Bonomo et al.| ( |201^aj ), in spite of the very poor sampling of the 
phase-folded LC transit (see their Fig. 5). The transit duration 
derived with SC data is slightly shorter (at 1.3 cr) than found by 
[Bonomo et al.[ ( [2012a[ ) but with negligible infiuence on system 
parameters. The planet eccentricity is consistent with zero within 
2 cr. 

From the stellar rotation period Frot = 12.95 + 0.25 d derived 
with all the LC data, the gyrochronology age t^y^ = 4 

agrees well with that estimated from stellar evolutiona ry tracks 
2.3 ^q^ Gyr, as already noted by Bonomo et al. (2012a). 


4.2.4. Kepler-44 

New system parameter s slightly differ from those derived by 
Bonomo et al. ( 2012a[ ) (see Table [^. Our analysis with SC 
data reveals a transit duration that is 2cr shorter than found by 
Bonomo et al.[ ( 2012^ . This implies a slightly larger < 2 //?★, hence 
a moderately higher stellar density from the Kepler third law. 
In consequence, stellar evolutionary tracks point to a slightly 
smaller star than previously reported by [Bonomo et al. ( [2012a[ ), 
with radius and mass of 1.35 + 0.08 R© and 1.12 + 0.08 M©, 
respectively. Therefore, also the planet turns out to be smaller 
(from Rp/R^) and denser: R^ = 1.09 + 0.07 Rju^ Mp = 1.00 + 
0.10 Mjup, andpp = 0.93 ^q|^ g cm“^ (see Tablej^. 

Figure shows the phase-folded transit and radial-velocity 
curve and, superimposed, the transit and the Keplerian models. 
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Fig. 3. Left panel, phase-folded transit light curve of Kepler-4 lb along with the transit model (red solid 
line). Right panel phase-folded radial-velocity curve of Kepler-41 and, superimposed, the Keplerian model 
(red solid line). 
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Fig. 5. Left panel phase-folded transit light curve of Kepler-43b along with the transit model (red solid 
line). Right panel phase-folded radial-velocity curve of Kepler-43 and, superimposed, the Keplerian model 
(red solid line). 


4.2.5. Kepler-74 

The circular solution we decided to adopt differs from the eccen¬ 
tric system parameters by almost three standard deviations. This 
is mainly because the transit density used to determine stellar 
parameters is a function of the eccentricity. For a null eccentric¬ 
ity, the planet becomes smaller and denser, with R^ = 0.97 + 
0.04 Rjup, Mp = 0.64 + 0.10 Mj^p, andPp = 0.86 + 0.18 g cm“^. 
The best fit of the transit and radial velocities is displayed in 
Fig.|7] 


The transit and radial velocity data along with the best solu¬ 
tions are shown in Fig.[^ 


4.2.7. KOI-205 


As for Kepler-75, the SC orbital and p hysical parameters agree 
very well with those determined by Diaz et al. ( |2013| ) (see 
Table [^. The best fit of the SC transit and radial-velocity ob¬ 
servations is shown in Fig.[^ 


4.2.6. Kepler-75 


The agreement between the system parameters determined by 
[Hebrard et al.| ( |2013'] ) and our DE-MCMC solution obtained with 
SC data (see Table is excellent, even adopting the slightly 
different atmospheric parameters derived with MOOG: Teff = 
5200 + 100 K and [Fe/H] = 0.30 + 0.12. 


The stellar rotation period inferred from the whole LC light 
curve is Fro t = 19.18 + 0.15 d, in agreement with Hebrard] 
|et al. (|2013|). The system age estimated from gyrochronology 
( [Mamajek & Hillenbrand|200^ is 1.6+0.3 Gyr, which is slightly 
lower than the value provided by stellar models, although con¬ 
sistent with the latter at 1.7 cr. 


5. Discussion and conciusions 

The analysis of Kepler SC photometry, RV data, SOPHIE and 
ESPaDONs spectra of seven giant companions has permitted us 
to refine their orbital and physical parameters. In three cases, 
namely Kepler-43, Kepler-75, and KOI-205, they agree with 
published parameters within 1 - 1.3 cr. 

For Kepler-44, for which only two quarters of LC data (Q1 
and Q2) were analysed for the discovery announcement, the new 
transit parameters determined with our DE-MCMC approach in¬ 
dicate a transit duration shorter by 2 cr, hence a slightly larger 
ajRi, and higher stellar density. This, in turn, implies that the 
host star and its planetary companion are smaller than previously 
found. 
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Fig. 6. Left panel, phase-folded transit light curve of Kepler-44b along with the transit model (red solid 
line). Right panel phase-folded radial-velocity curve of Kepler-44 and, superimposed, the Keplerian model 
(red solid line). 
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Fig. 7. Left panel phase-folded transit light curve of Kepler-74b along with the transit model for a circular 
orbit and their residuals. Right panel phase-folded radial-velocity curve of Kepler-74. Red diamonds and 
blue circles show SOPHIE and HARPS-N radial velocities, respectively. The black solid line displays the 
Keplerian circular model. 


A separate discussion must be made for Kepler-39b and 
Kepler-74b because we have revised the significance of their or¬ 
bital eccentricities. That of Kepler-39b is detected with a 2 cr 
significance level and, according to the Lucy-Sweeney criterion, 
it might be spurious. Slight asymmetries in the RV curve caused 
by residual effects from the correction of moonlight contami- 


nation (Santerne et al.| 

2011b 

) and/or CCD charge transfer in- 

efficiency (|Bouchy et a 

□2009 

1 might cause false eccentricities. 


Indeed, these effects become strong when observing faint stars. 
More radial-velocity observations without moonlight contami¬ 
nation are required to determine whether Kepler-39b has a low 
eccentricity. Kepler-74 would also benefit from additional high- 
accuracy and high-precision RVs because the chains of our DE- 
MCMC combined analysis did not converge towards a unique 
solution when including the eccentricity as a free parameter. Eor 
this reason, we decided to fit only a circular model to Kepler 
and RV data. Eor circular orbits, both Kepler-39b and Kepler- 
74b would be smaller and denser than previously found. In any 
case, the radius of Kepler-39b is still larger than predicted by 

theoretical isochrones (see Eig.Hl), as highlighted bylBouchy et| 

- 

Our new spectral analyses of the available co-added spectra 
revealed slightly hotter effective temperatures of Kepler-39 and 


KOI-205 and significantly higher metallicities for Kepler-39 and 
Kepler-75. 

Among our seven targets, the most striking divergence with 
already published parameters was found for Kepler-41. Indeed, 
SC data point to a considerably lower inclinat ion, higher im- 
pact parameter, and lower ajR^ than found by |Santerne et af 


( 2011a ) and |Southworth| ( [2()12| ). In consequence, both the host 
star and the planet have larger radii than previously derived. This 
new solution also has an impact on the estimation of the plane¬ 
tary geometric albedo that is significantly lower than previously 
estimated: Ag < 0.135 . Both the larger radius and the lower 
albedo make this planet resemble the majority of hot Jupiters. 
Conversely, the analysis of Kepler EC data had erroneously re¬ 
sulted in peculiar characteristics for this planet. This emphasizes 
that, in some cases, SC data are necessary to derive accurate 
system parameters , in addition to reducin g covariances between 
transit parameters ( Price & Rogers and permitting to com¬ 
pute precise transit timing variations. This will also be consid¬ 
ered in the light of the future TESS and PLATO space missions. 
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Fig. 8. Left panel, phase-folded transit light curve of Kepler-75b along with the transit model (red solid 
line). Right panel phase-folded radial-velocity curve of Kepler-75. Red diamonds and blue circles show 
SOPHIE and HARPS-N radial velocities, respectively. The black solid line displays the Keplerian model. 
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Fig. 9. Left panel phase-folded transit light curve of KOI-205 along with the transit model (red solid line). 
Right panel phase-folded radial-velocity curve of KOI-205 and, superimposed, the Keplerian model (red 
solid line). 
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Table 2. Kepler-39 system parameters. 


Fitted system parameters 

[Bouchy et al. (|2011 

) This work (eccentric) 

This work (circular) 

Orbital period P [days] 

21.0874 ± 0.0002 

21.087210 + 0.000037 

21.087212 + 0.000030 

Transit epoch ro[BJDTDB - 2454900] 

72.5959 + 0.0006 

1042.60708 + 0.00024 

1042.60707 + 0.00021 

Transit duration Tu [h] 

6.02 + 0.09 

5.960 + 0.020 

5.960 + 0.016 

Radius ratio Rb /R* 

Inclination i [deg] 

0.0896!«:“|‘ 

QQ QQ+0.59 

OO.OJ Q4Q 

U.UViU_QQQQg 

89.07 + 0.22 

0.0911 +0.0006 
89.23:°;}^ 

A 99+0.06 

033-8:?^ 

'^••^•^-0.09 

Limb-darkening coefficient qi 
Limb-darkening coefficient q 2 


n 072+0.08 
■^-0 06 

0 32+0-14 
-0.10 

sfe cosu; 

- 

-0 047+^-^^'^ 

0 (fixed) 

sfe sinu; 

- 

0 

u.jz.^-0 132 

0 (fixed) 

Orbital eccentricity e 

0.122 + 0.023 

0.112 + 0.057 

0(fixed) 

Argument of periastron [deg] cjj 

98.9!^:^ 

99+22 

^^-14 

90 (fixed) 

Radial velocity semi-amplitude K [kms"^] 

1.251 + 0.030 

1.257 + 0.064 

1.201 + 0.050 

Systemic velocity Vj- [kms"^] 

1 

p 

o 

1 + 

o o 

b b 

-0.063 + 0.044 

-0.032!«;“« 

RV jitter [ms“^] Sj 

- 

108!“ 

140!^« 

Derived transit parameters 




a/R^ 

Stellar density [g cm"^] 

23.8+;« 

u.j / _Q 11 

24.92!;| 

0 66+° 

-Q.ll 

27.74«« 

0 91+0.06 
^•^^-0.05 

Impact parameter b 


0 36+° °^ 

-0.08 

n 727+0.04 
-0.06 

Limb-darkening coefficient Ua 

0.303 + 0.014 « 

0.31 + 0.07 

0.32 + 0.06 

Limb-darkening coefficient 

0.308 + 0.005 « 

0.17 + 0.14 

0.15 + 0.11 

Atmospheric parameters of the star 




Effective temperature reff[K] 

6260 + 140 

6350 + 100 

6350 + 100 

Spectroscopic surface gravity log g [cgs] 

4.1 +0.2 

4.40 + 0.15 

4.40 + 0.15 

Derived surface gravity log g [cgs] 

4.19 + 0.07 

4.25 + 0.06 

4.34 + 0.02 

Metallicity [Fe/H] [dex] 

-0.29 + 0.10 

0.10 + 0.14 

0.10 + 0.14 

Stellar rotational velocity V sin u [km s"^] 

16 + 2.5 

16 + 2.5 

16 + 2.5 

Spectral type 

F8IV 

F7V 

F7V 

Stellar and planetary physical parameters 




Stellar mass [M©] 

-0.06 

f^9+0I}6 

^•^^-0.07 

1 26+°-°^ 

-0.06 

Stellar radius [R©] 

1.39 + 0.11 

1.40 + 0.10 

1.25 + 0.03 

Companion mass Mb [Mjup] 

18.0 + 0.9 

20 1+^-^ 

-Q.IO 

13.0+1° 

19.1 + 1.0 

Companion radius Rb [Rjup] 

Companion density pb [g cm"^] 

1 99+0.12 

12 4+^-'^ 

1.11 + 0.03 

17 4+^-^ 

4 /•+-1.4 

Companion surface gravity log gb [cgs] 

4.48 + 0.09 

4.51 + 0.05 

4.58 + 0.03 

Age t [Gyr] 

5.1 + 1.5 

9 1 +0.8 
^•^-0.9 

l-O-av 

Orbital semi-major axis a [au] 

0.155 + 0.003 

0.164 + 0.003 

0.162 + 0.003 

Equilibrium temperature Teq [K] ^ 

905 + 39 

897 + 29 

853 + 14 


^ The limb-darkening coefficients were allowed to vary within their Icr errors related to the uncertainties on stellar 
atmospheric parameters. 

^ Black-body equilibrium temperature assuming a null Bond albedo and uniform heat redistribution to the night side. 
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Bonomo et al. 2014: Improved parameters of seven Kepler giant companions 


Table 3. Kepler-41 system parameters. 


Fitted system parameters 

[Santerne et al. 

<2011af 

This work 

Orbital period P [days] 

1.S5555S + 0.00000/ 

1.85555820 + 0.00000052 

Transit epoch ro[BJDTDB - 2454900] 

70.1803 + 0.0003 

287.280482 + 0.000051 

Transit duration T^ [h] 

2.376 + 0.048 


2.4341 + 0.0046 

Radius ratio Rp/R^ 

0.0895 + 0.0019 

0.10265 + 0.00042 

Inclination i [deg] 

88.3 + 0.7 


82.51+0.09 

Limb-darkening coefficient qi 

- 


0.39 + 0.04 

Limb-darkening coefficient q 2 

- 


oo 

O O 

d d 
+ 1 

cn 

d 

Orbital eccentricity e 

0 (fixed) 


0 (fixed) 

Radial velocity semi-amplitude K [ms"^] 

85 + 11 


84+11 

Systemic velocity Vr [kms"^] 

-27.066 + 0.007 

-27.066 + 0.008 

RV jitter [ms"^] Sj 

- 


16!f 

Derived transit parameters 




a/R, 

6.43 + 0.05 


5.159 + 0.023 

Stellar density [g cm"^] 

1.46 + 0.04 


0.754 + 0.010 

Impact parameter b 

0.19 + 0.07 


0.672 + 0.005 

Limb-darkening coefficient Ua 

0.57 + 0.13 


0.46 + 0.07 

Limb-darkening coefficient 

0.41+0.13 


0.16 + 0.10 

Atmospheric parameters of the star 




Effective temperature reff[K] 

5620 + 140 


5750 + 100 

Spectroscopic surface gravity log g [cgs] 

4.20 + 0.15 


4.20 + 0.10 

Derived surface gravity log g [cgs] 

4.47 + 0.12 


4.278 + 0.005 

Metallicity [Fe/H] [dex] 

0.29 + 0.16 


0.38 + 0.11 

Stellar rotational velocity V sin U [ km s"^] 

6 + 2 


6 + 2 

Spectral type 

G6V 


G2V 

Stellar and planetary physical parameters 




Stellar mass [M©] 

1.12 + 0.07 


1.15 + 0.04 

Stellar radius [R©] 

1.02 + 0.03 


1.29 + 0.02 

Planetary mass Mp [Mjup] 

0.55 + 0.09 


0.56 + 0.08 

Planetary radius R^ [Rjup] 

0.89 + 0.05 


1.29 + 0.02 

Planetary density Pp [g cm"^] 

1.10 + 0.18 


0.33 + 0.04 

Planetary surface gravity log gp [cgs] 

3.23 + 0.09 


2.92 + 0.06 

Age t [Gyr] 

-0.3 


4 4+1-3 

Orbital semi-major axis a [au] 

0.030 + 0.010 


0.03101 + 0.0004 

Equilibrium temperature Teq [K] ^ 

1730 + 40 


1790 + 31 


^ Black-body equilibrium temperature assuming a null Bond albedo and uniform heat redistribution 
to the night side. 
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Bonomo et al. 2014: Improved parameters of seven Kepler giant companions 


Table 4. Kepler-43 system parameters. 


Fitted system parameters 

[Bonomo et al. 

(|2012af 

This work 

Orbital period P [days] 

5.024095 + 0.000021 

3.02409309 + 0.00000020 

Transit epoch ro[BJDTDB - 2454900] 

165.4159 + 0.0006 

195.45227 + 0.00005 

Transit duration T 14 [h] 

2.926 + 0.019 


2.900 + 0.003 

Radius ratio Rp/R* 

0.0868“ 

84.35!«:^J 


0.08647 + 0.00019 

Inclination i [deg] 


84.57!«:i? 

Limb-darkening coefficient qi 

- 


0.36 + 0.02 

Limb-darkening coefficient ^2 

- 


0.30 + 0.04 

sfe COS6U 

- 


0.061 + 0.030 

sfe sincj 

Orbital eccentricity e 

< 0.025 


0.086!0|J 

-0.009 

Argument of periastron [deg] oj 

- 


52+^^ 

^^-81 

Radial velocity semi-amplitude K [ms"^] 

375 + 13 


373 + 9 

Systemic velocity 14 [kms"^] 

-37.591 + 0.007 

-37.591 + 0.007 

RV jitter [ms"^] Sj 

- 


<6 

Derived transit parameters 




alR, 

Stellar density [g cm"^] 

Impact parameter b 

^01 +0.24 

-0.03 


6 93+^-ii 

0.648 + 0.004 

Limb-darkening coefficient Ua 

0.375 + 0.026 

a 

0.36 + 0.04 

Limb-darkening coefficient 

0.277 + 0.015 

a 

0.24 + 0.05 

Atmospheric parameters of the star 




Effective temperature reff[K] 

6041 + 143 


6050 + 100 

Spectroscopic surface gravity log g [cgs] 

4.64 + 0.13 


4.4 + 0.1 

Derived surface gravity log g [cgs] 

4.26 + 0.05 


4.26 + 0.02 

Metallicity [Fe/H] [dex] 

0.33 + 0.11 


0.40 + 0.10 

Stellar rotational velocity V sin U [ km s"^] 

5.5 + 1.5 


5.5 + 1.5 

Spectral type 

F 8 V 


F 8 V 

Stellar and planetary physical parameters 




Stellar mass [M©] 

1.32 + 0.09 


1.27 + 0.04 

Stellar radius [R©] 

1.42 + 0.07 


1 OO+0.05 

Planetary mass Mp [Mjup] 

3.23 + 0.19 


3.13 + 0.10 

Planetary radius [Rjup] 

1.20 + 0.06 


-0.03 

Planetary density Pp [g cm"^] 

2.33 + 0.36 


9 40+0.16 

'^-0.03 

Planetary surface gravity log gp [cgs] 

3.75 + 0.04 


Age t [Gyr] 

0 Q+1.0 


9 9+0.8 
^••^-0.7 

Orbital semi-major axis a [au] 

0.0449 + 0.0010 

0.0444 + 0.0005 

Equilibrium temperature T^q [K] ^ 

1637 + 47 


1628 + 33 


^ The limb-darkening coefficients were allowed to vary within their Icr errors related to the uncer¬ 
tainties on stellar atmospheric parameters. 

^ Black-body equilibrium temperature assuming a null Bond albedo and uniform heat redistribution 
to the night side. 
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Bonomo et al. 2014: Improved parameters of seven Kepler giant companions 


Table 5. Kepler-44 system parameters. 


Fitted system parameters 

[Bonomo et al.j(|2012a|) 

This work 

Orbital period P [days] 

5.246/40 + O.OOOOIS 

3.2467293 + 0.0000030 

Transit epoch ro[BJDTDB - 2454900] 

166.3781 + 0.0004 

287.15640 + 0.00021 

Transit duration T^ [h] 

3.218 + 0.043 

3.124 + 0.014 

Radius ratio Rp/R^ 

0.0844 + 0.0011 

0.0828 + 0.0008 

Inclination i [deg] 

Limb-darkening coefficient qi 

83.78!«;36^ 

84.96^ 

Limb-darkening coefficient q 2 

- 

0.33:0.15 

sfe cosu; 

- 

0 05+0-12 

sje sindu 

- 

0.03 + 0.20 

Orbital eccentricity e 

< 0.021 

< 0.066 

Radial velocity semi-amplitude K [ms"^] 

124 + 5 

125 + 11 

Systemic velocity 14 [kms"^] 

-35.892 + 0.004 

-35.892 + 0.009 

RV jitter [ms"^] s-^ 

- 

<7 

Derived transit parameters 



alR, 

Stellar density [g cm"^] 

6 45+^-^^ 

o:48|; 

' -0.04 

7 07+0.35 

'•^'-0.37 

0.63 + 0.10 

Impact parameter b 

0.62 + 0.02 

Limb-darkening coefficient Ua 

0.420 + 0.026 ^ 

0.42 + 0.13 

Limb-darkening coefficient Ub 

0.248 + 0.017 ^ 

0.21 + 0.19 

Atmospheric parameters of the star 



Effective temperature refF[K] 

5757 + 134 

5800 + 100 

Metallicity [Fe/H] [dex] 

0.26 + 0.10 

0.15 + 0.10 

Spectroscopic surface gravity log g [cgs] 

4.59 + 0.14 

4.1 + 0.1 

Derived surface gravity log g [cgs] 

4.15 + 0.06 

4.22 + 0.04 

Stellar rotational velocity V sin U [ km s"^] 

4 + 2 

4 + 2 

Spectral type 

G2IV 

G2IV 

Stellar and planetary physical parameters 



Stellar mass [M©] 

1.19 + 0.10 

1.12 + 0.08 

Stellar radius [R©] 

1.52 + 0.09 

1.35 + 0.08 

Planetary mass Mp [Mjup] 

1.02 + 0.07 

1.00 + 0.10 

Planetary radius Rp [Rjup] 

1.24 + 0.07 

1.09 + 0.07 

Planetary density Pp [g cm"^] 

0.65 + 0.12 

0 93+019 

Planetary surface gravity log gp [cgs] 

3.21 + 0.05 

3.31 + 0.06 

Age t [Gyr] 

6.95!}:^ 

5.8!?;^ 

Orbital semi-major axis a [au] 

0.0455 + 0.0013 

0.0446 + 0.0011 

Equilibrium temperature Teq [K] ^ 

1603 + 51 

1544 + 47 


^ The limb-darkening coefficients were allowed to vary within their Icr errors related to the uncer¬ 
tainties on stellar atmospheric parameters. 

^ Black-body equilibrium temperature assuming a null Bond albedo and uniform heat redistribution 
to the night side. 
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Bonomo et al. 2014: Improved parameters of seven Kepler giant companions 


Table 6. Kepler-74 system parameters. 


Fitted system parameters 

[Hebrard et al.j(|2013j) 

This work 

Orbital period P [days] 

7.540/18 + 0.000001 

7.340711 + 0.000006 

Transit epoch ro[BJDTDB - 2454900] 

67.3453 + 0.0003 

287.56737 + 0.00014 

Transit duration Tu [h] 

2.699 + 0.038R 

3.082 + 0.011 

Radius ratio R^/R^ 

0.090 + 0.002 

0.0912 + 0.0009 

Inclination i [deg] 

85.55 + 0.96 

87.46 + 0.07 

Limb-darkening coefficient qi 

- 


Limb-darkening coefficient q 2 

- 


Orbital eccentricity e 

0.287 + 0.062 

0 (fixed) 

Argument of periastron o) [deg] 

64 + 21 

- 

Radial velocity semi-amplitude K [ms“^] 

58 + 7 

59 + 11 

HARPS-N systemic velocity Vj-hn [kms"^] 

19.356 + 0.008 

19.366 + 0.009 

HARPS-N RV jitter [ m s"* ] ^j.hn 

<5 

12+16 

SOPHIE systemic velocity 14,so [kms"^] 

1 Q 9QO+0.008 

19.290 + 0.014 

SOPHIE RV jitter [ m s"^] 5j,so 

28 + 13 

40+15 

^^-10 

Derived transit parameters 



a/R^ 

11 .8+*'^ 

0.58!«;22 

15.47 + 0.18 

Stellar density [g cm"^] 

-0.04 

Impact parameter b 

0.684 ± 0.032 

0.685 + 0.011 

Limb-darkening coefficient Ua 

-0.17 

0.35 + 0.17 

Limb-darkening coefficient 

0.6 + 0.4 

^•+'^-0.23 

Atmospheric parameters of the star 



Effective temperature refF[K] 

6050+ 110 

6000 + 100 

Spectroscopic surface gravity log g [cgs] 

4.2 + 0.1 

4.5 + 0.10 

Derived surface gravity log g [cgs] 

4.2 + 0.1 

4.44 + 0.01 

Metallicity [Ee/H] [dex] 

0.34 + 0.14 

0.42 + 0.11 

Stellar rotational velocity V sin h [ km s"^] 

5.0+ 1.0 

5.0+ 1.0 

Spectral type 

E8V 

E8V 

Stellar and planetary physical parameters 



Stellar mass [M©] 

“T40!™ 

1.18 + 0.04 

Stellar radius [R©] 

1.51 +0.14 

1.12 + 0.04 

Planetary mass Mp [Mjup] 

0.68 + 0.09 

0.63 + 0.12 

Planetary radius R^ [Rjup] 

1.32 + 0.14 

0.96 + 0.02 

Planetary density Pp [g cm"^] 

0.37 + 0.13 

0.88 + 0.18 

Planetary surface gravity log gp [cgs] 

2.98 + 0.11 

3.23 + 0.10 

Age t [Gyr] 



Orbital semi-major axis a [au] 

0.084 + 0.014 

0.0781 + 0.0007 

Equilibrium temperature Teq [K] ^ 

1250 + 120 

1078 + 19 


This v alue of transi t duration is not accurate enough because it was estimated with the approximat ed formula 
reported in|Winn|(|2010|). Indeed, the transit duration is not a free parameter in the transit model adopted by|Hebrard et| 
^( [20T3] ). 


Black-body equilibrium temperature assuming a null Bond albedo and uniform heat redistribution to the night side. 
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Bonomo et al. 2014: Improved parameters of seven Kepler giant companions 


Table 7. Kepler-75 system parameters. 


Fitted system parameters 

[Hebrard et al. 

(|2013|) 

This work 

Orbital period P [days] 

S.SS4924 + 0.( 

J00002 

8.8849116 + 0.0000034 

Transit epoch ro[BJDTDB - 2454900] 

102.9910 + 0.0002 

840.44030 + 0.00012 

Transit duration T 14 [h] 

1.872 + 0.025 

u 

2.0778 + 0.0084 

Radius ratio R^jR^ 

0.121 + 0.002 


n 1 91 9+0.0009 
'^•^^^^-0.0007 

89.12!«;5> 

Inclination i [deg] 



Limb-darkening coefficient qi 

- 


0.36 + 0.08 

Limb-darkening coefficient ^2 

- 


0 42+010 
^•^^-0.08 

sfe cosu; 



0.337 + 0.020 

sfe sincj 



0.676 + 0.013 

Orbital eccentricity e 

0.569 + 0.010 


0.570 + 0.010 

Argument of periastron cj [deg] 

63.6 + 1.4 


63.5 + 1.7 

Radial velocity semi-amplitude K [kms"^] 

1.288 + 0.024 


1.287 + 0.025 

HARPS-N systemic velocity 14, hn [kms"^] 

-64.175 + 0.050 

-64.180 + 0.056 

HARPS-N RV jitter [ m s“^] ^j,HN 

< 17 


<54 

SOPHIE systemic velocity 14,so [kms"^] 

-64.235 + 0.012 

-64.236 + 0.014 

SOPHIE RV jitter [ m s“^] ^j,so 

< 8 


< 8 

Derived transit parameters 




a/R^ 

19.6 + 0.6 


1 Q 77+0.38 
' ^-0.45 

1 QC+O.Il 

L89_p 12 

Stellar density p„ [g cm"^] 

1.79 + 0.17 


Impact parameter b 

0.14 + 0.14 


0 13+0-09 

Limb-darkening coefficient Ua 

0.53 + 0.09 


0.50 + 0.05 

Limb-darkening coefficient 

0.13 + 0.26 


0.09 + 0.12 

Atmospheric parameters of the star 




Effective temperature reff[K] 

5330 + 120 


5200 + 100 

Spectroscopic surface gravity log g [cgs] 

4.55 + 0.14 


4.60 + 0.15 

Derived surface gravity log g [cgs] 

4.5 + 0.1 


4.50 + 0.02 

Metallicity [Ee/H] [dex] 

-0.07 + 0.15 


0.30 + 0.12 

Stellar rotational velocity V sin [ km s"^] 

3.5 + 1.5 


3.5 + 1.5 

Spectral type 

G 8 V 


KOV 

Stellar and planetary physical parameters 




Stellar mass [M©] 

0.88 + 0.06 


0.91 + 0.04 

Stellar radius [R©] 

0.88 + 0.04 


0.89 + 0.02 

Planetary mass Mp [Mjup] 

9.9 + 0.5 


10.1 + 0.4 

Planetary radius R^ [Rjup] 

1.03 + 0.06 


1.05 + 0.03 

Planetary density Pp [g cm"^] 

11+2 


ii.ot^ 

Planetary surface gravity log gp [cgs] 

4.36 + 0.03 


4.36 + 0.03 

Age t [Gyr] 

6 + 3 



Orbital semi-major axis a [au] 

0.080 + 0.005 


0.0818 + 0.0012 

Equilibrium temperature at the averaged distance Tgq [K] ^ 

850 + 40 


767 + 16 


^ This value of transit duration is not accurate enough because it was estimated with the approximated 


formula reported in 

Winn 

2010 

1 . Indeed, the transit duration is not a free parameter in the transit model 

adopted by 

Hebrarc 

et al. 

2013 

1 . 


^ Black-body equilibrium temperature assuming a null Bond albedo and uniform heat redistribution 
to the night side. 
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Bonomo et al. 2014: Improved parameters of seven Kepler giant companions 


Table 8. KOI-205 system parameters. 


Fitted system parameters 

[Diaz et al. <2013J 

This work 

Orbital period P [days] 

11.7201248 + 0.0000021 

11.720126 + 0.000011 

Transit epoch ro[BJDTDB - 2454900] 

75.17325 + 0.00012 

286.13609 + 0.00016 

Transit duration T^ [h] 

3.07 + 0.15 

3.069 + 0.017 

Radius ratio Rp/R^ 

0.09849 + 0.00049 

0.09906 + 0.00094 

Inclination i [deg] 

88.46 + 0.05 

88.43 + 0.06 

Limb-darkening coefficient qi 

- 

0 45+^-^^ 

-Q02 

Limb-darkening coefficient q 2 

- 

-0.15 

sfe cosu; 

- 

0.00 + 0.05 

^Je sindu 

- 

-0.03 + 0.010 

Orbital eccentricity e 

< 0.010 

< 0.015 

Radial velocity semi-amplitude K [kms"^] 

3.732 + 0.039 

o 7 C 7 + 0.071 
- -0.084 

Systemic velocity 14 [kms"^] 

15.057 + 0.026 

14.922 + 0.050 

RV jitter [ms"^] s-^ 

~40 

120 +^^ 

Derived transit parameters 



a/R, 

25.20 + 0.42 

25.07 + 0.43 

Stellar density [g cm"^] 

2.18 + 0.10 

2.17 + 0.11 

Impact parameter b 

0.676 + 0.014 

0.690 + 0.011 

Limb-darkening coefficient Ua 

0.523 + 0.017 ^ 

0.47 + 0.17 

Limb-darkening coefficient 

0.187 + 0.011 ^ 

0.21 + 0.24 

Atmospheric parameters of the star 



Effective temperature reff[K] 

5237 + 60 

5400 + 75 

Spectroscopic surface gravity log g [cgs] 

4.65 + 0.07 

4.7 + 0.1 

Derived surface gravity log g [cgs] 

4.550 + 0.015 

4.558 + 0.014 

Metallicity [Fe/H] [dex] 

0.14 + 0.12 

0.18 + 0.12 

Stellar rotational velocity V sin i^ [ km s"^] 

2 .0+ 1.0 

2.0 + 1.0 

Spectral type 

KOV 

G7V 

Stellar and planetary physical parameters 



Stellar mass [M©] 

0.92 + 0.03 

■096!^^ 

Stellar radius [R©] 

0.84 + 0.02 

0.87 ±0.02 

BD mass Mb [Mjup] 

39.9 + 1.0 

40.8!1:> 

BD radius Rb [Rjup] 

0.81 + 0.02 

0.82 ± 0.02 

BD density pb [g cm“^] 

100.4 + 6.9 

90.9!^:^ 

Age t [Gyr] 


i.rii 

Orbital semi-major axis a [au] 

0.0987 + 0.0013 

0.1010 + 0.0010 


^ The limb-darkening coefficients were allowed to vary within their Icr errors related to the uncer¬ 
tainties on stellar atmospheric parameters. 
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